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1 Introduction 



The experimental analysis of diffractive vector meson production by virtual or quasi- 
real photons (7J9 ^ Vp*) at HERA 00131111310111113131111! has been accompanied 
by numerous theoretical and phenomenological studies, e.g. [HI US [13 CH 113 113 113 
113 113 123 im 123 1211 • The main questions under discussion are the typical ones for 
semi-hard processes: To what extend perturbative QCD applies? What can we learn 
about non-perturbative hadronic interactions? 

The factorization proof gives a first answer. In the helicity amplitudes with 
^iil) = factorization holds and the contribution of a small-size qq dipole coupled 
by two gluons to the exchanged pomeron dominates by power counting the contribu- 
tions with additional soft exchanges. For the remaining helicity amplitudes the power 
counting does not result in the dominance of the short distance contribution. This is 
different from the related deep virtual Compton amplitude (DVCS, 7*p l*p) where 
the short distance contribution dominates in all helicity amplitudes (26j . 

In some presentations the situation is commented by calling the dominating short 
distance contribution in the helicity amplitude Aj = A/ = the leading twist one and 
by saying that short distance factorization must not be expected for the other helicity 
amplitudes because of being of non-leading twist. 

Owing to the DVCS case the latter argument looks not convincing. Moreover, the 
amplitudes with helicities different from A, = A/ = are accessible to experiment, e.g. 
in the ratio of the cross sections of longitudinal to transverse virtual photons and in the 
decay-angular distributions parametrized by the Schilling- Wolf ratios 123 C]- If short 
distance factorization holds the amplitude can be expressed in terms of generalized 
parton distributions of the target and predictions for ai/cTT [13 113 113 113 120] and 
the Schilling- Wolf ratios [23 I22j have been obtained. In some of the approaches the 
amplitudes have been treated in analogy to the one with Aj = Aj = 0, i.e. with 
the qq dipole interacting via two gluon exchange coupled to the generalized gluon 
distribution. However, unlike the Aj = A/ = case, in the leading twist contribution 
to the other helicity amplitudes singularities in the momentum fraction of the quark 
in the dipole appear. They are related to a large transverse size of the dipole and seem 
to signal the factorization breakdown expected from the power counting analysis |25j . 

There are controversial opinions about the appropriate treatment of those end- 
point singularities: introducing physically motivated cut-offs, choosing damping meson 
wave functions or including damping quark formfactors, or relating the argument of 
the parton distributions to the increasing dipole size. In any case, the partial success 
of these perturbative approaches seem to result in a modification of the first answer 
to the above main questions as given by the factorization proof: The factorization 
breaking effects cannot be large for all helicity amplitudes. 

The connection of the end-point singularities to the factorization breaking is re- 
ferred to frequently. In the study [27j using the operator product language these 
end-point singularities appear as the obstacle to factorization. 

In the present paper we end up with the opposite conclusion: The end-point con- 
tributions are factorizable; the ones appearing in the quark dipole interacting by two 
gluon exchange are factorized by identifying them as a leading \nQ^ Bjorken evolution 
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term of the two gluon {gg )to a quark - anti quark {qq ) exchange. Our discussion 
rehes on a specific model of the meson hght-cone wave function. It includes more than 
the leading twist contribution in terms of distribution amplitudes. Although finally 
only a small part of the specific information encoded in this wave function enters the 
results on the large asymptotics the inclusion in the wave function and resumma- 
tion of a geometric series of higher twist terms ~ (^if)'^ is the essential point for 
understanding the physical meaning of the end-point contributions. 

In section 2 we start with the impact factor representation of the diffractive ampli- 
tude and specify to the contribution of a scattering qq dipole coupled to the pomeron 
exchange by two gluons. In section 3 the '~f*V impact factors are introduced by speci- 
fying the meson light cone wave function. The asymptotics for large is calculated. 
We compare with the 7*7* impact factor and continue this comparison in the fol- 
lowing. Some details of the calculations are given in the Appendix. In section 4 
the logarithmic end-point contributions to the amplitudes are identified with a lead- 
ing InQ^ contribution to the generalized Bjorken evolution. In order to make this 
main point clearer this identification is repeated for the 7*p —>■ 7*p amplitude with 
Xi = Xf = 1, i.e. for the well known case related to the structure function F2{x,Q'^). 
In section 6 we evaluate numerically the resulting leading twist terms of the helicity 
amplitudes, specifying a the small x parton distribution, and obtain results for quan- 
tities which can be compared with the results of the experimental analysis: the cross 
section ratio aj^/aT and the angular-decay distribution (Schilling- Wolf) coefficients r"^ 
in dependence on Q"^ and t. 



2 Diffractive 7*V amplitudes 

A good starting point for analyzing high energy diffractive processes is the impact 
factor representation of the corresponding amplitude written in terms of partial waves, 

^"'"'M.i) = s r ^F^''{u;,Q,q) \(-—^ — + ( ) 

J-ioo27rz l\M\Q,m,q)J \M\Q,m,q) J 

f'''^ (to, Q,q) = j (fK^K'^'''^ {k, Q, q) g{K, k', q, u) ^^{k', q) (2.1) 

This is a typical form obtained in perturbative analysis PHUSHI; however it is based 
on more general arguments relying on impulse approximation jSU]. The field repre- 
senting the exchange interaction acts on the scattering particles for a short time, much 
shorter than the time scale of their binding or self-interaction. The field sees just a 
short-time intermediate fiuctuation state. The impact factors appear as the transition 
matrix elements (7* — >■ V and P — > P* in our case ) of an operator representing the 
action of the exchange field on the intermediate fiuctuations selecting some of the 
fiuctuations according to their interaction strength. 

The diffractive exchange Q is called Pomeron. In phenomenology it is often sub- 
stituted by a Regge pole or, according to a resent proposal [2H], by two poles, a soft 
and a hard Pomeron. In the framework of QCD we suppose the Pomeron to consist 
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out of interacting gluons. In perturbative QCD this idea acquires a definite meaning 
in the BFKL scheme [Ml ES]. 

In the case of hard diffraction the intermediate fluctuations is squeezed into a 
narrow space-time region (~ ^) in the vicinity of the light cone. Only short-distance 
modes of Pomeron field can interact with this fluctuation state. In this situation 
both the fluctuation state and the coupling to the Pomeron can be represented by 
perturbative QCD. Moreover one finds, that the fluctuations with a small number of 
partons dominates, the higher Fock states being suppressed by powers of ^. Also 
the exchange with a minimal number of exchange partons dominates not only by 
small couplings but, more important, by powers of ^. This suppression holds if the 
exchanged partons carry large transverse momenta k , ^ ^ m^, and one has 
to make sure that contributions with extra soft exchange quanta are absent or can be 
absorbed into generalized parton distributions. 

In the diffractive vector meson production by virtual photons the factorization in 
terms of a qq intermediate fluctuation coupling to the pomeron by two gluons has 
been proven on a rigourous level in the case where the short distance scale is provided 
by the momentum squared of the virtual photon with longitudinal polarization |,25,. 

In the case of factorization via two exchanged partons one can write 

gik, k', q, iu) = ^^^i— k', q, to) (2.2) 

In the considered case of electro-production we have in eq. ()2.1I ) M'^{Q'^,m,q) 
~ Q^. We pick up the leading contribution in the u integral by writing 



= J £K%{K,q,Q)j^-^^^^—^G'p{x^,X2,q, k) 



(2.3) 



G'p stands for the unintegrated generalized parton distribution of the proton (for 
recent reviews see |36[ ITTl ) resulting from the convolution of the proton impact 
factor with the projection of the exchange Q that couples by two exchange partons 
(p) in the small Bjorken variable limit, Xi = X2 = The skewdness ^ = ^J-=^ 
is small; we have non-vanishing transverse momentum transfer q. 

As a simplification we replace by the derivative of the gluon distribution at 
small X 

G; = ^j^[G,(xi,X2,|g|M«:r) T{k',Q'))] (2.4) 

Here T(k^, Q"^) is the parton Sudakov formfactor; it is equal to 1 at = and small 
for < Q^. 

The relation to the standard notation introduced by Ji is as follows, 

rJi/ .^ X1 + X2 ^ Xi- X2 



G,{xr,X2) = Hi\x,0, x=-^^,i=-^^^. (2.5) 
For the leading contribution in we expand $p'^^(fc, g, Q"^) 

^'''^ (., g, q) = (^^^^) ^ ^ '^^^2'*^' ^ c^^'^ + ... 



(2.6) 
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The leading In contribution of the k integral is then obtained as 



(2.7) 



As an obstacle to factorization in diffractive electroproduction with transverse 
polarization, mentioned in the Introduction, one encounters large contributions from 
qq fluctuation states with the longitudinal momentum fraction of one of the partons 
{z or z = 1 — z) small. The soft scattering quark then may couple by soft exchange 
gluon to the Pomeron without suppression by 1/Q. 

We shall investigate the expressions for impact factors with qq state and a model 
wave function ip^ of the vector meson. We shall show that the enhanced end-point 

I 1 2 

contribution at large {s ^ ^ my) actually arise from 1 ^ z,z -f^. In 
this range the z integral is approximately logarithmic and this contribution can be 
identified as the one of generalized Bjorken evolution [221 122] of the two exchanged 
gluons gg to exchanged quark-antiquark qq. The latter exchange involve higher twist 
modes. 

In this way we are going to show that the soft quark or end point contribution 
can be factorized and included in the parton distribution of the diffractive exchange. 
The factorization of this contributions goes via a qq exchange instead of gg. The gg 
coupling to the Pomeron is then restricted to the contribution of the hard (2;, z = 0(1) 
) scattering qq dipole. 

3 Impact factors 

The impact factor with qq intermediate state and two leading gluon exchange can be 
written as (compare Fig. 1) 



The first argument in the light-cone wave functions ipi/f is the transverse momen- 
tum relative to the momentum direction of corresponding particle. are the trans- 
verse momenta components of the exchange gluons, Ki + = —q is the transverse 
part of the momentum transfer. 

The 7*7* impact factor can be treated purely perturbative. It is obtained in 
the perturbative Regge asymptotics as the integral over the right hand cut of the 
discontinuity of the amplitude 7* (7 7* (7. The four terms in (j3.1|) correspond 

to the four ways to couple the two gluons to gg; one of them is shown in Fig. lb. The 
momentum variables are defined in the Sudakov frame 




2) = y dH^dH^dz^jl^ih, ^)</)^*^'°'^(£i, /2, Ki, K2W/*{l2 - Zq, Z) 

0"*P°'^(£i, £2, ^^2) = Ois[5\l2 - ^1) + 5\(i2 - ^1 + «:i + /^a)- 

5\^2-^l + ^^l)-5\^2-^l + ^^2)] (3.1) 



Qi = q' - xip, Q2 = q' ~ X2P + g, 
l = zq' - I3f^p + £, 
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Figure 1: a) Impact factor form of the 7*P VP amplitude. The rectangular box in- 
dicates the exchange interaction and the dashed line box the unintegrated generalized 
gluon distribution, b) Contribution to the 7*F impact factor. 
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ki = aiQ - (3ip + Ki, 
2q'p = s 



(3.2) 



The 6 functions of the mass shell condition at the right-hand discontinuity in the 
subenergy —2kiQi ^ Pis, e.g. for Fig. lb, 

z6{{z - ai){pi - pe)s - |£ - - mj) z5{z{(3i - xi)s - - mj) (3.3) 

are used to do the integrals over this subenergy (3is and also over the loop momentum 
Sudakov component Pe. The Sudakov components ai of the exchanged gluons are 
neglected in this step. 

The two remaining propagators result in the light-cone wave function of 



7' 



'^^'^^ (£, z, Q) e |/i|2 1 ™2 

[Q' + 

^(0) ^ yi+1) ^ ^(-1) ^ ^ 

z z 



(3.4) 



The sum over fermion chiralities adds the same with z ^ z or amounts in the factor 
2 in dSIH) 

In the helicity cases A, = A/ = ±1 there is an additional chiral odd term which is 
included by the substitution 



1 1 TTl 

2v^\e„ z)v*^\h, z) = tM-, + + 

As a model for the light cone wave function of the vector meson we assume 



^^^(£,z)=/, 



V^{i,z,mv) 



m 



exp 



V 



zzm 



V J 



(3.5) 



(3.6) 



The form is motivated by QCD sum rules, it is formally obtained by Borel transfor- 
mation of the propagator factor in ()3.4|) with respect to and by the substitution 
of Borel variable by my, where my is of the oder of the meson mass [3T]. This wave 
function, being close to the one of 7*, is a particular realization of the phenomeno- 
logically successful concept of vector dominance. Actually the explicit form of 
involves more information than necessary for the asymptotic estimate at large Q. The 
essential point in changing from to is removing the hard (singular in impact 
parameter, the Fourier conjugate to €) component while keeping the helicity structure. 

In fl3.4l 13. 6p we have kept the quark mass to indicate the possible extension to the 
case of heavy quark vector mesons; it will be neglected in the following. 

In the Appendix we consider the impact factors in some detail and calculate the 
leading twist contribution, for 7*7* at Q\ = xis, Q\ = X2S, with xi,X2 small and 
s — > cx) and for 7*V^ at — 00, for representative cases of polarizations. For the 
asymptotic estimate we divide the range in z into zq < z < 1 — zq = zq, < z < zq 
and Zq < z < 1,, 



(3.7) 
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The results for z = 0(1) are 



l-zo 



<^^'^^ (g, k)= ifl'^^' (g, k)zzdz, 



zz^Tiq,k)=Cr^^-^+ 



I /■(2,**)|2 1 

fraclzz - (4 - -) 
2 ^ Q4 ^ ^ zz^ 



z-z^fiq,k) = C''i^-^iA-^), (3i 

zz 

zM'i<l,k)=C^'^^^-^{A-^), 



z-z^\'\q.k) = Cl''l^^ + 

Cl'"^Ay'\k,q)f'>*{k,q)^+[q*'f'\k,q)~\q\'f^'**)^ (4-3)1, 

.^^-(g,fc) = 0-^^^(1-2). 

We have introduced the abbreviations 

/(2)(A;, g) = k{k + g)* + k*{k + g), /(2)**(fc, g) = 2k*{k + g)*, 

/(=')(fc,g)=g7('nfc,?) + ^g/^'^**- (3.9) 

In the case 7*1^ the large scale is just Q^. For a smoother extrapolation back 
into the subasymptotic region we replace = + my. In the case 7*7* the large 
scale is to be substituted in ()3.8|) by = s; actually s enters always multiplied 
by Xiy + which we have absorbed into the coefficients C^^^f in this case. 

The coefficients C^^^^ depend on xi, X2 in the case of 7*7* : 

C7°'°(xi X2) = -2^ X1X2Y {xi, X2, 1,1,2), C2°'°(xi X2) = 2^/x^Y{xl,X2, 2, 2, 3) 

C^°''(xi,a;2) = v^r(xi,X2,2,l,2) 

C^i'0(xi X2) = -V2Y{xi, X2, 1, 2, 2) (3.10) 
Cr''\x,X2) = 2Y{x,, X2, 1, 1, 1), Cr''\x,X2) = -2Y{x^,X2, 2, 2, 2) 
C^i'i(xiX2) =X2r(xi,a;2,0,2,2), 

with the abbreviation 



Y{xi, X2, a, 6, c) 



{xiy + X2y) ' 



In the case 7*V^ the coefficients f depend on mv and Q 

^1 — *-'2 — ^T^T o~ 



C 



V0,1 _ f^vQ 



V 

2 



+ g2)3/2 



^ -(m2. + g2)i/2 ^^-^^^ 



-<yi,-i 1^^1,-1 r^Vl.l 



In the asymptotically dominating amplitude Aj = = we have included the 
next-to leading term in the twist power expansion. This is done in oder to demonstrate 
that the higher twist terms are accessible in the present approach and to test their 
effect on the numerical estimates. 

The leading terms of the end-point contributions from z = O{zo) are 

_ c^'^^ in ^ - in ^ + 0(.„), (3,12) 

g" \f^? g" \q? 

The contributions from the other end point, ^zo, are obtained by replacing zq by 
Zq. Here Cq'^^ coincide with the coefficients in front of J ^ in in eq. (|3.8|) as 

expected for cancelling the auxiliary zq. 



4 Factorization of end point contributions 

As the result of the previous section we have separated the end point contributions. 
Starting from the impact factor with a reasonable meson wave function we have iden- 
tified the end point contributions for the kinematics ^ ~ |gp as being 
proportional to 




We observe that there is no divergence at z, z^ 0. Spurious divergences appear only 
in the extrapolation of the twist expansion done for z = 0(1) (13.81) to z,z^ 0. As 

2 

we see, the blind extrapolation would neglect terms with powers of -^^qi- The model 
vector meson wave function just specifies how the sum of these higher twist terms 
removes the end point divergence and leads to ( 14.11) . This point would be missed 
when looking only at the leading twist term of the qq dipole scattering with two gluon 
coupling to the exchange. 

The 7*7* impact factor shows the same structure of end point contributions at 
large s. The amplitude is the one of non-forward virtual Compton scattering in the 
Regge asymptotics xi,X2 <C 1. The standard short distance factorization for the 
Compton amplitude at xi,X2 = 0(1) can be continued into this region. In leading 
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In we have two cases. The hard scattering sub-process may be the one of Compton 
scattering off a quark which enters the amphtude in convolution with the (generahzed) 
quark distribution. This means, the coefficient functions starts at tree level; the quark 
loop contributes only in next-to-leading order. In the other case the hard scattering 
sub-process is the one of 7*5' '~f*g via a quark loop entering the amplitude in 
convolution with the (generalized) gluon distribution; the corresponding coefficient 
function starts with the one-loop order. This suggests the way how to treat the end 
point contributions in the 'j*V case. 

We start from the limiting form impact factor at large from which we have 
obtained the end-point contributions ()3.12|1 (compare also Appendix ()8.5p8.12j) ). 

$^0 = ^ / e ^'"v / dz 

mf, J z Jo 



Q"^ + qn* + q*K + — -|- 



...iK + q = 0)... } (4.2) 



We have restored the integration over the loop transverse momentum i ()8.4|1 . i' = 
i — Y^i^, + zq) ^ i — K. Further we restore the integrations over the loop momentum 
Sudakov components Pi, (|3.2|) . and over the subenergy (gi — fci)^ ^ PiS by including 
the mass shell 5-functions fl3.3|) . 

= / sdf3i / d(3e ... z6{z{(3i - /3^)s - |f p) 6{{Pe - xi)s - \i' + (4.3) 



^0 



The periods stand for the right-hand side of ()4.2|1 . Now the first 5-function is used 
to do the integral over z; in particular we do the substitution = {(3i — f3e)s. We 
notice that the first term in the argument of the second (5-function dominates for large 
Q"^ = xis. 

We retain only the contributions with a logarithmic contribition of the z integral. 
If substituted into the amplitude only the structure f^'^^n, q) results in a logarithmic k 
integral allowing to approximate the proton impact factor by the unintegrated gener- 
alized gluon distribution. In fact we pick up the leading contribution in the kinematics 

2 

V- 
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2 



Jo Jo 

2e '"V 

The resulting end-point contribution to the amplitude has the form 

i f(")fr/i r^^ d^f r^^ r^^ 




(4.4) 



Jml 



Jo Jo J\k\'^ 

P^(Ai,A2;/?i,/52) G,(A,/?2,g;K'r), (4.5) 
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where 



This is just the evolution term in the integral form of the generalized GLAPD equation, 
corresponding to the parton splitting gg — > qq. The splitting kernel, 

mi,/5£2;/5i,/32) = iPi-Pii)-' exp(-^i^^) [1 + ^1^]-!-, (4.6) 
is an unconventional one, because the resulting qq state is not the one of leading twist 

2 

exchange. We have X2 = I3(,2 = — ^ <^ 1 and therefore the kernel results in a narrow 
distribution peaked around = Pi. 



5 Generalized Bjorken evolution 

We consider in some details the end-point contribution of the 7*p '~f*p amplitude 
with Xi = Xf = 1. We expect to recover the GLAPD evolution gg qq known from 
the case of the structure function F2 and the non- forward generalization (DVCS). 

For simplicity we restrict ourselves to purely longitudinal momentum transfer, 
g = 0, Xi 7^ X2. The end-point contribution to the impact factor reads 



<l>^'i'i(A;,0)| 



d^i I dz 




* 1 



XiS + ^ 



X2S + 



\e-K\ 



^ -...(«: = 0)... 



(5.1: 



The subtraction term is obtained from the written term by substituting k = 



^71 in HM 2 / dH' f\ ndz 

^^'^'\k,0)\,^, = e'as I / dy — 

^ Jo Jo z 



\n -y\f^\ 



z 



^0 



s{xiy + X2y) + + H-yy 



2 I d^^' f\ r dz \ -2^-^\K\^yy-y\K\h{xiy + X2y) + 0{\ti\^^ 
e as / / dy 



z 



s{xiy + X2y) + 



As in fj3.3p we restore the integration over f3i, f3e by including the mass shell S function. 
This allows to substitute in the integrand 

If P 

= - /?ijs = (/?2 - ^2)*, Ai = a;i,A2 = ^2 



and in particular 



s{xiy + X2y) + 



l\2 



We obtain in analogy to the previous section 



$^'i'i(A;,0)|,„ = e^ 



s{Piy + p2y) 



d(3id(32d(3i^5{xi - X2 - Pi + (32 
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S{Ph-xi) dy 0/?i-Xi 5.3 

Jo [/3iy + l32y\^ 

The result contributes to the amphtude 7*p 7*p (Aj = A/ = 1) by convolution with 
the unintegrated gluon distribution, 



v 



GgiP^,P2;q = 0,\i\') (5.4) 

As expected, the result has the form of the evolution term in the integral representation 
of the GLAPD equation with the (non-forward) splitting kernel 

P{x^X2■, PuP2) = + 0{f3, - X,)) (5.5) 

P1P2 

We compare P{xiX2; Pi, P2) with the generalized GLAPD evolution kernel as calcu- 
lated in [40j 

H{xi,X2,Pl,p2) = ^{a;i(^l - ^11') + (xi -/?i)Ji + (/?i ^/?2)} 

00 

Jl = [axi + 1 — i£]^^[—ax2 + 1 — i£:]~""^[a(xi — /?i) + 1 — ie]^^ 

J -00 27r2 ^ 

Jii' = / [axi + 1 — ie]~^[a{xi — f3i) + 1 — ie]~^ (5.6) 

The first (second) term in the bracket corresponds to the gg — > qq transition with- 
out (with) spin flip. For the impact factor the right-hand cut discontinuity in Pis is 
relevant 

rr 1 (XiX2 {Xl - (3l){X2 - p2)\ n ^ 

= Aft Ift ft + ftft - 

= P(l,,l2;ftft) (5,7) 

The familiar forward splitting kernel for gg — >■ qq is recovered for jSi = l32,xi = X2- 

The comparison of (5.7) and (5.5) confirms that the end-point contribution to the 
diffractive amplitude written in impact factor representation is to be absorbed into 
the evolution of the generalized gluon distribution. Since the impact factor ansatz 
accounts for the Regge asymptotics the corresponding evolution contribution results 
merely in this approximation; in particular the subenergy {Qi — kiY — s{f3i — Xi) is 
small in this asymptotics. 



6 Amplitudes and numerical estimates 

The decomposition of the z range in the impact factor ()3.7|) results in the amplitude 
as a sum of corresponding 3 terms, M'^^'^f = M'^"^* + M^J''^-' + M^^'^* . According to 
fl2.3p the contribution from 2; = (9(1) is calculated as 

A A Z"*^^ Z"^" A A 1 

M^" = / d^K I zz(p^" ^{q,K)dz ^^(xi, X2, g, |ft:p). (6.8) 

Jmlr Jzo \^\ \q + '^l 
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The leading In contribution of the integral over k results in 

mI"^' ^ (^j^ ° zzq^l'^'U^ ■ X2, g; Q'). (6.9) 

Here ^^^^ are the coefficients in (j3.8|) accompanying f^'^\n-,q) as in (|2.6p . For 
example, in the case Aj = A/ = 1 we have 

/ (_-2) rfz , ■Gg(xi,X2,g;g^). (6.10) 

J zo / V + 

The end point contributions are not small with zq for the terms in ()3.8|1 involving ^. 
As explained in Section 4 the logarithmic z integral is a generalized Bjorken evolution 
()4.5p term resulting in the effective quark distribution 

G,(a;i,X2,g,zoQ')= T"^ ^'^^^G,{x,,x,,qM''). (6.11) 
In the example Aj = A/ = 1 this leads to 

= = n2^^2 ■ G',(xi,X2,g,^oQ') (6.12) 

The sum of the three contributions to the amplitude has the form (j2.7p . Notice 
that the effective quark distribution has been indentified with the Bjorken evolution 
term of the gluon distribution. This involves the reasonable assumption that this 
distribution is small for small x at the scale m^. 

The generalized gluon distribution at small Xi/2 is adopted here to be proportional 
to the ordinary small x gluon distribution 

G,(xi, X2, g, Q') = c Gg(^^^, Q') e""'^'' (6.13) 

For the slope parameter we adopt the value 6 = 6 GeV"^. 

The generalized Bjorken evolution leads to a modification of this relation if it is 
assumed to hold at some Q^. This effect has been investigated for small x and small 
skewedness ISH] and has been applied to diffractive vector meson production In 
the latter paper the dependence of the cross sections is found to be changed by 
50 per cent for light vector mesons at the highest (= SOGeV"^). Since the aim 
of the numerical estimates in this paper is merely the illustration of the proposed 
factorization concept, we do not include the skewdness effects here for simplicity, 
understanding that this would be an appropriate further improvement. 

Further, for the numerical estimate we have to specify the gluon distribution. Al- 
though our approach is based on the asymptotic expansion in for comparison with 
the data we have to extrapolate the results to non-asymptotic values of this scale. 
Moreover, in the reconstruction of the effective quark distribution according to 1)6.111) 
we have to integrate the small x gluon distribution starting from my ~ .5GeV^. This 
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means we need the gluon distribution at small scales, where the standard parametriza- 
tions like MRST do not give a certain answer and where the application of the evolution 
equation is not reliable. Therefore we adopt the two-pomeron paramerization 



2 n2 



where Eq = 0.43, ei = 0.08, Xq = 0.0014, Xi = 0.5954, Ql = 9.108, Qf = 0.5894. 

Again, we have no particular reason to favour this parametrization besides of its 
convenience in our aim of illustration. We did not try to optimize the results with 
regard to different possible input parametrizations. 

The cut-off zq is to separate the z range of order unity from the end-point regions; 
a value zq = 0.1. ..0.2 is reasonable and indeed in this range the zq dependence of the 
result is weak. The condition that zqQ'^ is much larger than my leads to the restriction 
of the applicability of the estimates to > 5GeV^. 

We have calculated in the given approximation with the input specified above the 

dependence of the ratio of longitudinal to transverse virtual photon diffractive 
cross sections. We do this including all helicity amplitudes (R) and also with the 
helicity conserving amplitudes only (Rq), the latter case corresponds to the ratio 
usually obtained in the data analysis. Thus we calculate 

where 

a^(Q2) = j dt{\M''\^ + |Mi°|2 + iM^-^n, 4°)(Q2) = j dt\W\^. (6.15) 

The results are shown in Fig. 2 together with HERA data. Rq deviates from R by 
10 per cent at higher Q^. The inclusion of the next-to- leading terms in the twist 
expansion does not lead to noticeable changes. 

The coefficients in the angular-decay distribution are more sensitive to the helicity 
dependence, because in some of them the small flip amplitudes enter in the first 
power. We use the following expressions for the coefficients r^^^ in terms of the helicity 
amplitudes M^^^f simplified in comparison to |21] for the appropriate case of the 
virtual photon polarization parameters being e ~ 1, 5 ~ 0. 
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Figure 2: Ratio of the longitudinal and transverse elastic p° electroproduction cross 
sections as a function of . The dotted line corresponds to the assumption of helicity 
conservation, the solid line takes into account spin flip amplitudes. 

N = |M°Y + + 2\M^'\^ + \M''\' + \M'-'\^ 

Results are given for Tqq, Tqq + 2rj'-|^, Tqq + 2rl^. In the data analysis the first 
coefficient is extracted from the dependence on the polar angle of the tt"*" in the vector 
meson rest frame with respect to the vector meson momentum direction. Actually Rq 
is obtained from this coefficient by a formula valid for the ratio in the approximation 
of helicity conservation. The other combinations arc the ones obtained in the data 
analysis from the dependence on the angle between the lepton scattering and vector 
meson production planes. 

The t dependence of the coefficients rfj has been calculated at a flxed value of 

= lOGeV'^. The results are shown in Fig. 3 in comparison with the data for a 
broad range = 3. ..30 GeV'^ with the average values of 6.5 GeV'^ for the ZEUS data 
points and 5 GeV'^ for the HI data points. 

In evaluting the dependence we have substituted in the amplitudes the average 
value of t determined by the slope parameter b. The results are shown in Fig. 4. 



7 Discussion 

The diffractive '-f*V amplitudes constructed with a qq dipole impact factor involve en- 
hanced contributions from the end-point regions in the quark longitudinal momentum 
fraction z, besides of the leading contribution to the Aj = A/ = helicity amplitude. 
These endpoint contributions appear as singularities in the extrapolation of the twist 
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Figure 3: Too + 2rfi, rgo + 2r^i, Tqq as a function of transferred momentum t. 



16 





;ure 4: Tqq + 2rfi, rl^ + 2r^i, r^Q as a function of 
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expansion . However, resummed corrections proportional to (^^)"' remove the sin- 
gularities at 2;, ^ = and result in logarithmically enhanced end-point contributions. 
The latter are shown to be the ones of the generalized Bjorken evolution of the t- 
channel parton exchange gg — > qq. The short-distance factorization of the diffractive 
amplitude in terms of the '~f*V transition and the Pomeron exchange therefore in- 
volves contributions of two types: one with a scattering qq dipole coupled by gg to 
the Pomeron and one with a (Compton hke) scattering g or g coupled by qq to the 
Pomeron. The qq exchange differs from the leading twist one. 

We have chosen the impact factor form of the diffractive ampltude as a starting 
point. Our lack of understanding of the proton impact factor and of the Pomeron 
coupling to it has been managed by replacing their convolution by the unintegrated 
generalized parton distribution, a standard step in kx factorization schemes. This is 
justified in so far as it results at large in the parton distribution convoluted with the 
coefficient function resulting from the 'j*V impact factor. Notice that in the leading 
In approximation contributions of the impact factors depending on the azimuthal 
angle of the exchanged parton momentum, -j^, drop out. We see here a possible source 
of corrections, which may be relevant at moderate Q^. 

The resulting large asymptotics of the amplitudes consists of terms with the 
generalized small x gluon distribution and with the effective qq distribution. The latter 
appears as an additional non-perturbative input. However, its main contribution re- 
sults from the evolution of the gluon distribution with a splitting function transferring 
the longitudinal momentum fraction from g to q almost unchanged. 

In this way we see that the construction of the diffractive amplitude with a scale 
of the gluon distribution replaced by zzQ"^ |Tn| is actually an approximation to the 
factorization proposed here, because with that scale replacement the end-point con- 
tributions of the ^-integration are approximately the qq exchange contribution by 
gg —>■ qq evolution. This scale replacement has been applied in a previous study of 
the polarization effects |21j . 

2 

In our construction the terms ~ (^qz)" improving the end-point behaviour are 
introduced via the vector meson light cone wave function. It is modelled starting from 
the perturbative 7* wave function by keeping its helicity structure and removing its 
hard contributions. Wave functions sharing these basic features lead to similar results 
for the large helicity amplitudes. In the present scheme the higher twist correction 
to the amplitudes can be calculated. The leading asymptotics involves the value of 
the impact parameter wave function at the origin and this can be recast in of the 
distribution amplitude formulation. 

We did not include the hard scale evolution of the wave function In the nu- 
merical estimates we did not include the skewedness effects. We have used a particular 
parametrization of the small x gluon distribution covering the small range needed 
in particular for reconstructing the effective qq distribution therefrom. 

The results for the cross section ratio and the Schilling- Wolf coefficients have been 
calculated using this input with no attempts of fitting. We have shown that they 
describe the basic features of the data. 

Including the mentioned improvements the proposed factorization scheme provides 
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a theoretical framework applicable to diffractive electroproduction free of previous 
uncertainties about the end-point contributions and therefore useful for extracting 
information about the structure of the hadrons involved from data analysis. We 
expect that the scheme applies to other hard diffractive processes. 
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8 Appendix 

We take the simple form of 1)3.111 into account and write 

$(/c,g)=/ dzzz[(p{k,z,q) + ip{-k - q,z,q) - (p{0,z,q) - ^{-q,z,q)], (8.1) 

^0 
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For the 7*7* case we have 



\i\r\r* 



\£—K—zq\^+m^ 1 



Doing the integration over the transverse momentum we obtain 



For the 'j*V case we have 



yy 



y y 



zzyy ~^ zz 



(8.3) 



exp(- 



1-1 
rriyzz 



The result of the integration over £2 can be written as 



ip^{K,z,q) = efvasTT 



1 + A 



1.5) 



Q\ m2(l + A) \K + zq\'' X 
exp[ — 



mr 



zzrnir 



rriyZZ 1 + A 



We hst the numerator {Vi^'Vj^^ (y, Qi, Q2) for the representative cases of hehc- 
ities: 



A,; 





1 
1 
1 



A/ 




-1 

-1 
1 



2Q1Q2 

-Q,y[n + zq){\-\) 
Q2y{n + zqY{\-\) 
yy{K + zqyj= 
-[\K + zq\^y + zzQl]{^ + ^) 



As expected, the resulting expressions obey 

(p{K,,z,q) = (p{-K - q,z,q) 
and in the integral (j8.ip the number of terms can be reduced to two. 



(8.6) 
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We calculate the asymptotics of the integral at first for z = 0{1) for the 7*7* case 
in the region s — > 00, = sxi, with fixed <^ 1 



s s'^ {xiy + X2yY zz s^' 

and for the case 7*1^ in the region ^ 00 



1.7) 



ip'^{K,z,q)^efvasTT{V/''Vf 
1 ln + zql"^ 



+ 0{- 



1 



Notice that the numerator includes terms proportional to s or if Aj = A/. A further 
term in the expansion has to be included, leading to cancellation if Aj = A/ = ±1. In 

2 

the case 7*1^ it is appropriate to choose the expansion parameter as q2^^2 resulting 
in an improvement of the extrapolation in the region of moderate > 
We obtain in the case 7*7* 



my. 



^'''\K,z,q)z-z = + Cr'ix„x,)^-^^±P^^ + 

s zz s-^ 

7,01/ \ - r^7,oi/ \k + zq\'^{k + zq) 1 1 ^/ 1 n 
(f^' {K,z,q)zz = C{' {xi,X2) ^ (--=) + 0(— ) 

7,10/ N - r^7,iO/ ,\K + zq\'^{K + zq)* l 1 /nr ^ \ 
{K,z,q)zz = CI' Xi,X2 {---) + 0{—) 

z, q)z-z = Cr'~\x^,X2) ^''^"^^*\ 

s 

c-'-'(x,..)<^i±£CJ^i±i«!!4 + o(A) 



^.9) 



zz 



(/^■^'^^(k, z, q)zz = C7'"(a;i, Xa) 



\k + zq\^ _ 1 



0(- 



S Z^ Z S"" 

We use the abbreviation for fl3.10|) . In the case 7*\^ we have similar fomulae. 
When calculating the integrand of (j8.1|) . 



V94 = V5(ft;, 2:, g) + V5(-K - g, ^, g) - V5(0, g) - ^{-q, z, g), 



^.10) 



the dependence on z and k, /t + g, disentangles. We list the relevant expressions in the 
following table 



^{K,z,q) 


V94(/€,2;,g) 


\k + zq\^ 
{k + zq)\K + zgp 

(K + ^g)*^ 
(k + zg)*^ K + 2;gp 
ft + -zg^ 


+ g)* + + g) = g) 
{z -z){qK{K + g)* + gK*(K + g)) = {z - z)f^^\K,q) 
2K*{K + qy = f^'**\K,q) 

l/(2)/(2,«) + (1 _ 3;,^)[g*2j(2) ^ |g|2j(2,**)] 
1|/(2,**1|2 + (1 _ 4^^)[|^|2J(2) ^ 1^*2J(2,**)* ^ lg2j(2,**)] 
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In this way we obtain the twist expansion vahd for z = 0(1). The result for 
the case 7*^ is given in eq. ()3.8p : the one for 7*7* can be recovered by substituting 

We calculate now the end point contribution to the asymptotic expansion in s for 

7*7* case, 

z'zLp'^^K, z,q)dz = dy dz 



[s{xiy + X2y) + yy{i^q* + i^*q + + -l^^LiE 
and in Q"^ for the 7*V case 



+ O{zo) (8.11) 



zzLp [K,z,q)dz= / dz 







[Q'^ + mlr + [nq* + K*q+\K\^) +^-^] 

+Oizo) (8.12) 



The numerator results from the expansion of z'z < VjVf >' for small z. Wi is non- 
vanishing only in the case Aj = A/ = ±1. In the case Aj = A/ = both Wq and Wi 
vanish, here the small z expansion in the numerator starts with W-i ~ z. 
Thus we are lead to calculate the integrals 

Ia{K,Q')= / dz (8.13) 
J, [Q2 + ^] 





\k\ 


2 






2 



' Q 

and obtain 



zzip'' (k, z, q)dz = Wi{y, k)Ii{k, Q') + 1^0(2/, «)/o(«;, Q ) + O{zo) 



In the case 7*7* we have to substitute by s{xiy + X2y) + yyi^q* + i^*q)', in the 

2 - — ^ 

vector- meson case we substitute y by Q2^i and by Q'^ + my + K*q + nq* . The 
small term {n*q + Kg*) matters in the case Aj = A/ = ±1 only. For illustration it is 
enough to do one case, 7*^, Aj = 0, A/ = 1, explicitely: Wi = 0, Wq = —yQn, 



|2 

In ^^^^ — "i—TTT — ' — "^'1. — In — ' — 7--^ ^'J!' In 
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At + q\ 
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K + g 
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Q4 |^ + g|2 Q, 

In the kinematics ZqQ"^ ^ ^ |gp this can be approximated by 
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2 



(m2 + Q2)3/2 [ (m^ + (^2)3/2 |^|2 j 

In this way we obtain ()3.12|) confirming the matching with the expression ()3.8p for 
z = 0{1), i.e. the cancellation of the auxiliary zq. 
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